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Abstract. The possibility to probe new physics scenarios of light Majorana neutrino exchange and right- 
handed currents at the planned next generation neutrinoless double f3 decay experiment SuperNEMO is 
discussed. Its ability to study different isotopes and track the outgoing electrons provides the means to 
discriminate different underlying mechanisms for the neutrinoless double /3 decay by measuring the decay 
half-life and the electron angular and energy distributions. 
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and that flavour is violated in the leptonic sector [T]. De- 
spite this success, oscillation experiments are unable to de- 
termine the absolute magnitude of neutrino masses. Upper 
limits on the effective electron neutrino mass of 2.3 eV [21 
and 2.05 eV [3| can be set from the analysis of tritium j3 
decay experiments. Astronomical observations combined 
with cosmological considerations yield an upper bound to 
be set on the sum of the three neutrino masses of the or- 
der of 0.7 eV [3]. However, the most sensitive probe of 
the absolute mass scale of Majorana neutrinos is neutri- 
noless double /3 decay (0v/3f3) [S[[SJ[71[S] . In this process, 
an atomic nucleus with Z protons decays into a nucleus 
with Z + 2 protons and the same mass number A under 
the emission of two electrons, 



(A, Z) -> (A, Z + 2) + 2e- 



(1) 



This process can be described by the exchange of a light 
neutrino connecting two V-A weak interactions, see Fig. [1] 
(a). The process ([1]) is lepton number violating and, in 
the standard picture of light neutrino exchange, it is only 
possible if the neutrino is identical to its own anti-particle, 
i.e. if neutrinos are Majorana particles. Combined with 
the fact that neutrino masses are more than five orders of 
magnitude smaller than the masses of other fermions, such 
a possibility suggests that the origin of neutrino masses is 
different from that of charged fermions. 

Several mechanisms of mass generation have been sug- 
gested in the literature, the most prominent example be- 
ing the seesaw mechanism [S] in which heavy right-handed 
neutrinos mix with the left-handed neutrinos and generate 
light Majorana masses for the observed active neutrinos. 
The Majorana character of the active neutrinos can then 
be connected to a breaking of lepton number symmetry 
close to the GUT scale and might be related to the baryon 
asymmetry of the Universe through the baryogenesis via 
leptogenesis mechanism |10j . 

Because of its sensitivity to the nature and magnitude 
of the neutrino mass, 0v/3(3 decay is a crucial experimen- 
tal probe for physics beyond the Standard Model and its 
discovery will be of the utmost importance. It will prove 
lepton number to be broken, and in most models it will 
also provide direct evidence that the light active neutrinos 
are Majorana particleqj [12) . However, the measurement 
of 0^/3/3 decay in a single isotope is not sufficient to prove 
that the standard mechanism of light Majorana neutrino 
exchange is the dominant source for the decay. There are a 
host of other models, such as Left- Right symmetry [5], R- 
parity violating Supersymmetry (SUSY) [T3J or Extra Di- 
mensions [11 , which can provide alternative mechanisms 
to trigger Oz//3/3 decay. In some of these models, additional 
sources of lepton number violation can supplement light 
neutrino exchange. For example, in Left-Right symmet- 
ric models, there are additional contributions from right- 
handed currents and the exchange of heavy neutrinos. In 
other models, such as R-parity violating SUSY, Oz//3/3 de- 
cay can be mediated by other heavy particles that are not 
directly related to neutrinos. 



See [TT] for a counter-example of a model where such a 
conclusion is not valid. 



There are several methods proposed in the literature 
to disentangle the many possible contributions or at least 
to determine the class of models that give rise to the domi- 
nant mechanism for 0^/3/3 decay. Results from Oz^/3/3 decay 
can be compared with other neutrino experiments and ob- 
servations such as tritium decay to determine if they are 
consistent. At the LHC there could also be signs of new 
physics exhibiting lepton number violation that is related 
to Oz^/3/3 (see [T3] for such an example in R-parity violat- 
ing Supersymmetry). Such analyses would compare results 
for 0^/3/3 with other experimental searches, but there are 
also ways to gain more information within the realm of 
Ov/30 decay and related nuclear processes. Possible tech- 
niques include the analysis of angular and energy correla- 
tions between the electrons emitted in the Oz//3/3 decay [SJ 
[T5l[T^l[T71[T5] or a comparison of results for Oi//3/3 decay in 
two or more isotopes [TW2TJll2"Tll2"2"] . These two approaches 
are studied in this paper. Other proposed methods are the 
comparative analysis of Oi//3/3 decay to the ground state 
with either Oz^/3 + /3 + or electron capture decay [231 and 
Oi//3/3 decay to excited states [24] . 

Currently, the best limit on Oi//3/3 decay comes from 
the search for 0^/3/3 decay of the isotope 76 Ge giving a 
half-life of T 1/2 > 1.9 • 10 25 years [25]. This results in 
an upper bound on the 0vf3(3 Majorana neutrino mass of 
(m u ) < 300 — 600 meV. A controversial claim of obser- 
vation of 0^/3/3 decay in 76 Ge gives a half-life of Ti/ 2 = 
(0.8-18.3) • 10 25 y [26 and a resulting effective Majorana 
neutrino mass of (m„) = 110 — 560 meV. Next genera- 
tion experiments such as CUORE, EXO, GERDA, MA- 
JORANA or SuperNEMO aim to increase the half-life ex- 
clusion limit by one order of magnitude and confirm or 
exclude the claimed observation. The planned experiment 
SuperNEMO allows the measurement of 0^/3/3 decay in 
several isotopes ( 82 Se, 150 Nd and 48 Ca are currently be- 
ing considered) to the ground and excited states, and is 
able to track the trajectories of the emitted electrons and 
determine their individual energies. In this respect, the 
SuperNEMO experiment has a unique potential to disen- 
tangle the possible mechanisms for 0^/3/3 decay. 

This paper addresses the question of how measure- 
ments at SuperNEMO can be used to gain information 
on the underlying physics mechanism of the 0^/3/3 decay. 
The sensitivity of SuperNEMO to new physics parame- 
ters in two models is determined by performing a detailed 
simulation of the SuperNEMO experimental set-up. By 
analysing both the angular and energy distributions in 
the standard mass mechanism and in a model incorporat- 
ing right-handed currents, the prospects of discriminating 
0^/3/3 decay mechanisms are examined. The two models 
are specifically chosen to represent all possible mecha- 
nisms, as they maximally deviate from each other in their 
angular and energy distributions. 

This paper is organised as follows. In Section [5] a short 
description of the theoretical framework on which the cal- 
culations of the Oz^/3/3 decay rate and the angular and en- 
ergy correlations are based is shown. The example physics 
models are introduced and reviewed. Section [3] gives a 
brief overview of the SuperNEMO experiment design and 
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Fig. 1: Diagrams illustrating 0vf3(3 decay through (a) the mass mechanism and (b) the right-handed current contribu- 
tion via the A parameter. 



in Section 0] a detailed account of the simulation anal- 
ysis and its results are presented. In Section [5j the ex- 
pected constraints from SuperNEMO on new physics are 
shown and the prospects of disentangling Oi>/30 mecha- 
nisms by analysing the angular and energy distributions 
and by comparing rates in different isotopes are addressed. 
Our conclusions are presented in Section [51 



2 Neutrinoless Double Beta Decay 

2.1 Effective Description 

Contributions to Qvj3(3 decay can be categorised as either 
long-range or short-range interactions. In the first case, 
the corresponding diagram involves two vertices which are 
both point-like at the Fermi scale, and connected by the 
exchange of a light neutrino. Such long-range interactions 
are described by an effective Lagrangian [2"Tll2"g] 




LA. 



lJv-A + 22 e abJaJb 



(2) 



a . b 



where Gf is the Fermi coupling constant and the leptonic 
and hadronic Lorentz currents arc defined as j a = eO a v 
and J a — uO a d, respectively. Here, O a denotes the corre- 
sponding transition operator, with a — V — A, V + A, S — 
P,S + P,T L ,T R [57]. In Equation ©, the contribution 
from V — A currents originating from standard weak cou- 
plings has been separated off and the summation runs over 
all Lorentz invariant and non-vanishing combinations of 
the leptonic and hadronic currents, except for the case 
a = b = V — A. The effective coupling strengths for long- 
range contributions are denoted as e l ^ b . 

For short-ranged contributions, the interactions are 
represented by a single vertex which is point-like at the 
Fermi scale, and they are described by the Lagrangian [251 

on 



r 2 
C = —f-m, 



L.i. 

1 V^ 

p / j *~abc 

a,b,c 



'J a'JbJ c' 



(3) 



Here, m p denotes the proton mass and the leptonic and 
hadronic currents are given by J a = uO a d and j' a = 
eO a e c , respectively. The transition operators O a are de- 
fined as in the long-range case above, and the summation 
runs over all Lorentz invariant and non-vanishing combi- 
nations of the hadronic and leptonic currents. The effec- 
tive coupling strengths for the short-range contributions 
are denoted as e*£ c . 

Described by the first term in Equation @, the ex- 
change of light left-handed Majorana neutrinos leads to 
the 0^/3/3 decay rate 



PI^T 1 = {{m v )/m e ) 2 G m \M mu f 



(4) 



where (mj) is the effective Majorana neutrino mass in 
which the contributions of the individual neutrino masses 
rrii are weighted by the squared neutrino mixing matrix 
elements, V\, (m u ) = | £, U^rm]. 

Analogously, other new physics (NP) contributions, of 
both long- and short-range nature, can in general be ex- 
pressed as 



K,?} 



e 2 NP G NP \M 



NP\ 



(5) 



where e^p denotes the corresponding effective coupling 
strength, i.e. is either given by e l ^ b for a long-range mech- 
anism or by e*£ c for a short-range mechanism. In Equa- 
tions (j5]) and (0 , the nuclear matrix elements for the mass 
mechanism and alternative new physics contributions are 
given by A4 m „ and SAnp, respectively, and Gqi, Gnp de- 
note the phase space integrals of the corresponding nuclear 
processes. It is assumed that one mechanism dominates 
the double (3 decay rate. 



2.2 Left-Right Symmetry 

The focus in this paper is on a subset of the Left-Right 
symmetric model [5., which incorporates left-handed and 
right-handed currents under the exchange of light and 
heavy neutrinos. Left-Right symmetric models generally 
predict new gauge bosons of the extra right-handed SU(2) 
gauge symmetry as well as heavy right-handed neutrinos 
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Isotope 


Is mm [y 


Caa [y" 1 ] 


C mX [y" 1 ] 


7b Ge 
82 Se 
150 Nd 


1.12 x 10"" 
4.33 x 10" 13 
7.74 x 10" 12 


1.36 x 10"" 
1.01 x 10" 12 
2.68 x 10" 11 


-4.11 x 10" 14 
-1.60 x 10" 13 
-3.57 x 10" 12 



matrix, and the light and heavy neutrino mass scales, m v 
and Mr, leads to the approximate relations: 



Table 1: Coefficients used in calculating the OV/3/3 decay 
rate [3D]. 



which give rise to light observable neutrinos via the seesaw 
mechanism. 

The OV/3/3 decay half-life in the Left-Right symmetric 
model can be written as a function of the effective param- 
eters /i,?7, A [50] , 



+C m \[i\ + C mr) \xr\ 



CxxX 2 + C m if 



(0) 



where contributions from the exchange of heavy neutrinos 
are omitted. The coefficients C mm ,C nn etc. are combina- 
tions of phase space factors and nuclear matrix elements. 
The first three terms give the contributions from the fol- 
lowing processes: 

1- C mm /Lt 2 : Fully left-handed current neutrino exchange, 
see Fig.[T](a) (mass mechanism); 

2. C\\\ 2 : Right-handed leptonic and right-handed had- 
ronic current neutrino exchange, see Fig. [T](b); 

3. C m ri 2 : Right-handed leptonic and left-handed hadronic 
current neutrino exchange. 

The remaining terms in Equation ([6]) describe interfer- 
ence effects between these three processes. The effective 
parameters /i, rj, A in ((SJ) are given in terms of the under- 
lying physics parameters as 



i=l 
3 

r^tanC^C/i, 1 ^, 



m r 



\Mw R 



2 

ei ' 



(7) 
(8) 
(9) 



i=l 



with the electron mass m e , the left- and right-handed W 
boson masses Mw L and Mw R , respectively, and the mix- 
ing angle £ between the W bosons. The 3x3 matrices 
U 11 and U 12 connect the weak eigenstates (i/ e , v^ v T ) of 
the light neutrinos with the mass eigenstates of the light 
neutrinos (v\, 1/2, V3), and heavy neutrinos, (Aq, N2, N3), 
respectively. We assume that the neutrino sector consists 
of three light neutrino states, m Vi <C m e , and three heavy 
neutrino states, Mpj i 3> m p , i = 1, 2, 3. Consequently, the 
summations in (O [HI El) are only over the light neutrino 
states. For a simple estimate of the sensitivity of OV/3/3 de- 
cay to the model parameters, we neglect the flavour struc- 
ture in U 11 and U 12 ; using the assumption that the ele- 
ments in U 11 are of order unity (almost unitary mixing), 
and those in U 12 are of order mo /Mr ~ y/m v /MR, with 
the effective magnitude mo of the neutrino Dirac mass 



/»' 


m v 
m e ' 




(10) 


n ~ 


itanC /S' 




(11) 


Ap 


. (Mw L \ 2 

" \m Wr ) V 


Mr 


(12) 



In the following analysis a simplified model incorporating 
only an admixture of mass mechanism (MM) due to a neu- 
trino mass term /x = (m v )/m e and right-handed current 
due to the A term (RHCa) is considered: 



P1/2] — C mm [i + C\\\ + C m \jj,\. 



(13) 



As we will see in Section l2Tl these two mechanisms exhibit 
maximally different angular and energy distributions, and 
with an admixture between them, to a good approxima- 
tion any possible angular and energy distribution can be 
produced. In our numerical calculation we use the values 
as given in Table [T] for the coefficients C mm , Caa and C m \ 
in Equation (|13p . Furthermore, we assume that the pa- 
rameter [i is real-valued positive and A is real- valued. 



2.3 Nuclear Matrix Elements 

As demonstrated in Equations ((3]) and ([5]), a calculation of 
the nuclear matrix elements (NMEs) is required to convert 
the measured half-life rates or limits into new physics pa- 
rameters. Exact solutions for the NMEs do not exist, and 
approximations have to be used. Calculations using the 
nuclear shell model exist for lighter nuclei such as 76 Ge 
and 82 Se, though the only reliable results are for 48 Ca. 
Quasi-particle random phase approximation (QRPA) cal- 
culations are applied for most isotopes as a greater num- 
ber of intermediate states can be included. In this paper, 
a comparison between two possible SuperNEMO isotopes 
( 82 Se and 150 Nd) and the isotope that gives the current 
best limit ( 76 Ge) is made. Consistent calculations of the 
NMEs for these three isotopes in both the MM and RHC 
are rare (only [30] and [31]). All the results are shown 
using NMEs from (3D], displayed in Tabled] 

Recent work on the calculation of NMEs for the heavy 
isotope 150 Nd suggests that nuclear deformation must be 
included, as QRPA calculations usually consider the nu- 
cleus to be spherical. To compensate for this a suppres- 
sion factor of 2.7 is introduced into the NME due to an 
approximation arising from the BCS overlap factor |32) . 
7W( 150 Nd)/2.7. This gives a suppression C mm ,AA,mA/(2.7) 2 
in Table [T] The 82 Se nuclei are assumed to be spherical 
and no correction is added in this paper. 

The NMEs are a significant source of uncertainty in 
double /3 decay physics and quantitative results in this 
paper could change with different calculations (particu- 
larly for 150 Nd). For example, more recent studies [33] 
suggest the NMEs from 150 Nd for the MM are an addi- 
tional factor 1.3-1.7 smaller. In our analysis we assume a 
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theoretical uncertainty of 30% in the NMEs of all isotopes 
and mechanisms considered throughout. Even though the 
choice of NME changes quantitative results for the ex- 
tracted physics parameters, the conclusions about the ad- 
vantages of using different kinematic variables will not be 
affected. 



2.4 Angular and Energy Distributions in the Left-Right 
Symmetric Model 

For our event simulation, the three-dimensional distribu- 
tion of the Of/3/3 decay rate in terms of the kinetic energies 
<i,2 of the two emitted electrons and the cosine of the angle 
between the electrons cos 6*12 is used: 



p(h,t 2 , cos 6»i 2 ) = 



dr 



dtidt2d cos 9i2 



(14) 



The distributions for the MM and for the RHCa mecha- 
nism are given by 

PMM(ti,t2,COS0i2) — 

d x (tx + l) Pl (t 2 + l)p 2 F(t 1 ,Z)F(t 2 , Z) 

x 5(Q - ti-t 2 ){l-Pip 2 cos 9 X2 ), (15) 

/9RHC(*1,*2,COS01 2 ) = 

c 2 x (ti + l)pi(t 2 + l)p 2 F{t 1 ,Z)F{t 2 , Z)(h - t 2 f 
xS(Q - t 1 -t 2 )(l + /3i/3 2 cos 12 ), (16) 



with the electron momenta pi — yti(ti + 2) and velocities 
Pi = Pi/(U + 1), an d the mass difference Q between the 
mother and daughter nucleus. All energies and momenta 
are expressed in units of the electron mass and c\ and 
c 2 are normalisation constants. The Fermi function F is 
given by 



F(t, Z)=c 3 x p 2s - 2 e™\r(s + iu)\ 2 



(17) 



where s = y/1 - (aZ) 2 , u = aZ(t + l)/p, a = 1/137.036, 
r is the Gamma function and C3 is a normalisation con- 
stant. Here, Z is the atomic number of the daughter nu- 
cleus. The normalisation of the distributions is irrelevant 
when discussing energy and angular correlations. 

Using Equations fli"5t and (fl"6|) . the differential decay 
widths with respect to the cosine of the angle 6\ 2 , 



dr 



-=/ dhpit^Q -h, cos 12 ), (18) 

12 Jo 



dcos6* 
and the energy difference At = t\ —t 2 

dr 



d{At) 



1 a a (Q±M Q-At 

d COS O12P , , COS #12 



(19) 
may be determined. 

The differential width in Equation (|18p can be written 
as [UTS] 

l^a =^(l-fc9cos0i 2 ), (20) 

d COS #12 2 



with the total decay width r. The distribution shape is 
linear in cos#i 2 , with the slope determined by the param- 
eter kg which can range between —1 < kg < 1, depending 
on the underlying decay mechanism. Assuming the domi- 
nance of one scenario, the shape does not depend on the 
precise values of new physic parameters (mass scales, cou- 
pling constants) but is a model specific signature of the 
mechanism. For the MM and RHCa mechanisms, the the- 
oretically predicted kg is found from Equation (fTS"]) and is 
given by 



k, 



So 






Se 



= +0.88, 
= -0.79, 



C'm = +0-89, 



(21) 
(22) 



The correlation coefficient kg is modified when taking into 
account nuclear physics effects and exhibits only a small 
dependence on the type of nucleus. The MM and the 
RHCa mechanisms give the maximally and minimally pos- 
sible values for the angular correlation coefficient kg in a 
given isotope, respectively. 

Experimentally, kg can be determined via the forward- 
backward asymmetry of the decay distribution, 



Ag = 



dr 

dcos9 



N+-N- 

N+ + N. 



■dcosO - 

kg_ 
2 ' 



dr 

dcost 



■dcose) IF 



(23) 



Here, N + (iV_) counts the number of signal events with 
the angle 612 larger (smaller) than 90°. 

Analogously, the MM and RHCa mechanism also differ 
in the shape of the electron energy difference distribution, 
Equation (TK)|) . For the isotopes 82 Se and 150 Nd, these dis- 
tributions are shown in Fig. [5] Again, the shape is largely 
independent of the isotope under inspection. The following 
asymmetry in the electron energy-difference distribution 
is determined, 




d(At) - 



Q/2 d(At) J 



(24) 



thereby defining an energy correlation coefficient &E, where 
Q is the energy release of the decay. The rate N + (iV_) 
counts the number of signal events with an electron en- 
ergy difference smaller (larger) than Q/2. For the MM and 
RHCa mechanism, the theoretical kE parameter may be 
found from Equation (|19l) and is given by 



4mm = +0-66, 



k 



So 



-1.07, 



■Nd 



<,Nd 



-0.64, 
-1.09, 



(25) 
(26) 



in the isotopes 82 Se and 150 Nd. As can be seen in Fig. [2j 
the MM and RHCa mechanisms correspond to different 
shapes of the energy difference distribution. Analogous to 
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Fig. 2: Normalised Oi/ftft decay distribution with respect to 
the electron energy difference in the MM (red) and RHCa 
mechanism (blue) for the isotopes 82 Se (solid curves) and 
150 Nd (dashed curves). 



the angular distribution, the corresponding energy corre- 
lation coefficients in the two mechanisms considered are, 
to a good approximation, at their upper and lower limits 
in a given isotope. 



3 SuperNEMO 

SuperNEMO is a next generation experiment building on 
technology used by the currently running NEMO-III ex- 
periment pi[551l551l571[551[5giHCJ] . The design of the de- 
tector consists of 20 modules each containing approxi- 
mately 5 kg of enriched and purified double ft emitting 
isotope in the form of a thin foil (with a surface density 
of 40 mg/cm 2 ). Isotopes under consideration for Super- 
NEMO are 82 Se, 150 Nd and 48 Ca. 

The foil is surrounded by a tracking chamber compris- 
ing nine planes of drift cells (44 mm diameter) on each 
side operating in Geiger mode in a magnetic field of 25 
Gauss. The tracking chamber has overall dimensions of 4 
m height (parallel to the drift cells), 5 m length and 1 m 
width (perpendicular to the foil) ; the foil is centred in this 
volume with dimensions of 3 m height and 4.5 m length. 
The tracking allows particle identification (e~,e + ,7,a) 
and vertex reconstruction to improve background rejec- 
tion and to allow measurement of double ft decay angular 
correlations. 

Calorimetry consisting of 25 x 25 cm 2 square blocks of 5 
cm thickness scintillating material connected to low activ- 
ity photomultiplier tubes (PMTs) surrounds the detector 
on four sides. An energy resolution of 7% (FWHM) and 



time resolution of 250 ps (Gaussian a) at 1 MeV for the 
blocks is required. The granularity of the calorimetry al- 
lows the energy of individual particles to be measured. Ad- 
ditional 7- veto calorimetry to identify photons from back- 
ground events of thickness 10 cm surrounds the detector 
on all sides. The modules are contained in shared back- 
ground shielding and will be housed in an underground 
laboratory to reduce the cosmic ray flux. A diagram of the 
planned SuperNEMO module design is shown in Fig. [3] 



4 Simulation 

4.1 Simulation Description 

A full simulation of the SuperNEMO detector was per- 
formed including realistic digitisation, tracking and event 
selection. Signals for two mechanisms of Qvftft decay (mass 
mechanism MM and right-handed current via the A pa- 
rameter RHCa) and the principal internal backgrounds 
were generated using DECAY0 [4"T]. This models the full 
event kinematics, including angular and energy distribu- 
tions. 

A GEANT-4 Monte Carlo simulation of the detector 
was constructed. Digitisation of the hits in cells was ob- 
tained by assuming a Geiger hit model validated with 
NEMO-III with a transverse resolution of 0.6 mm and a 
longitudinal resolution of 0.3 cm. The calorimeter response 
was simulated assuming a Gaussian energy resolution of 
7%/y/E (FWHM) and timing resolution of 250 ps (Gaus- 
sian a at 1 MeV) . Inactive material in front of the 7- veto 
was partially simulated. 

Full tracking was developed consisting of pattern recog- 
nition and helical track fitting. The pattern recognition 
uses a cellular automaton to select adjacent hits in the 
tracking layers. Helical tracks are fitted to the particles. 
Tracks are extrapolated into the foil to find an appropriate 
event origin and into the calorimeter where they may be 
associated with calorimeter energy deposits. The realistic 
event selection (validated using NEMO-III) was optimised 
for double ft decay electrons (two electrons with a common 
vertex in the foil). The selection criteria are: 

— events must include only two negatively charged par- 
ticles each associated with one calorimeter hit; 

— event vertices must be within the foil and the tracks 
must have a common vertex of <10 standard devia- 
tions between intersection points in the plane of the 
source foil; 

— the time of flight of the electrons in the detector must 
be consistent with the hypothesis of the electrons orig- 
inating in the source foil; 

— the number of Geiger drift cell hits unassociated with 
a track must be less than 3; 

— the energy deposited in individual calorimeter blocks 
must be > 50 keV; 

— there are zero calorimeter hits not associated with a 
track; 

— tracks must have hits in at least one of the first three 
and one of the last three planes of Geiger drift cells; 
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Fig. 3: A SuperNEMO module. The source foil (not shown) sits in the centre of a tracking volume consisting of drift 
cells operating in Geiger mode. This is surrounded by calorimetry consisting of scintillator blocks connected to PMTs 
(grey). The support frame is shown in red. 



— the number of delayed Geiger drift cell hits due to a 
particles from 214 Bi- 214 Po events must be zero; 

— there are no hits in the 7-veto detectors with energy 
> 50 keV. 

Using these experimental selection criteria the signal 
efficiency was found to be 28.2% for the MM and 17.0% 
for the RHC A in 82 Se and 29.1% for the MM and 17.3% 
for the RHCa in 150 Nd. This is higher than the efficiency 
for MM detection in 100 Mo decays in NEMO-III of 17.4% 
(in the electron energy sum range 2-3.2 MeV) |4T)] . 

The variables reconstructed from the simulation are 
the energy sum, where a peak is expected at the energy 
release, Q, of the Ov/3/3 decay, the energy difference and 
the cosine of the opening angle of the two electrons. Sim- 
ulations of the angular and energy difference distributions 
of the two electrons in a signal sample are shown in Fig. @] 
for the isotope 82 Se (similar results hold for 150 Nd). The 
reconstructed distributions, normalised to the theoretical 
distributions, show detector effects which arise due to mul- 
tiple scattering in the source foil, compared to the theo- 
retically predicted distributions based on Equations (fi"8l) 
and (|19p. This influence is particularly strong in the right- 
handed current as one electron usually has low energy so 
the shape of the distribution is changed (on average a 
30° deviation from the generated distribution). The recon- 
struction efficiency is also low for small angular separation 
between the electrons when they travel through the same 
drift cells. 

The backgrounds were processed by the same detector 
simulation and reconstruction programs as the signal. The 
dominant two neutrino double /3 decay (2z^/3/3) background 
and the background due to foil contamination were nor- 
malised assuming a detector exposure of 500 kg y. Due to 
the high decay energy Q for 0^/3/3 in 150 Nd, the 214 Bi back- 
ground may be neglected. The activities were assumed to 
be 2 /iBq/kg for 208 T1 and 10 /iBq/kg for 214 Bi. These 
are the target radioactive background levels in the base- 
line SuperNEMO design. Reconstructed distributions of 



the experimental variables including background events 
for the MM at an example signal half-life of 10 25 y are 
shown in Figure [5j 



4.2 Limit Setting 

To determine the longest half-life that can be probed with 
SuperNEMO, exclusion limits at 90% CL on the half- 
life using the distribution of the sum of electron ener- 
gies (Fig. [5] (a)) were set using a Modified Frequentist 
(CL S ) [42 method. This method uses a log-likelihood ra- 
tio (LLR) of the signal-plus-background hypothesis and 
the background-only hypothesis, where the signal is due 
to the Oz^/3/3 process. The effect of varying the 214 Bi back- 
ground activities on the expected limit to the MM is shown 
in Fig. El The expected limit is given by the median of 
the distribution of the LLR and the widths of the bands 
shown represent one and two standard deviations of the 
LLR distributions for a given 214 Bi activity. For com- 
parison, the NEMO-III internal 214 Bi background is < 
100 /iBq/kg in 100 Mo and 530 ± 180 ^Bq/kg in 82 Se. The 
NEMO-III internal 208 T1 background is 110 ± 10 /iBq/kg 
in 100 Mo, 340 ±50 ^Bq/kg in 82 Se and 9320 ±320 /iBq/kg 



l.iO 



Nd 



j . The 7-veto used reduces the number of ra- 



dioactive background events by 30% for 214 Bi in the elec- 
tron energy sum window > 2.7 MeV. 

All external backgrounds from outside the foil, apart 
from radon in the tracking volume, are expected to be 
negligible and were not simulated. The energy distribu- 
tion of the external radon background is similar to the 
internal background. Simulations have shown that a con- 
tamination of 10 /iBq/kg of 214 Bi in the foil is equivalent 
to 280 /iBq/m 3 of 214 Bi in the gas volume and 2 /iBq/kg 
of 208 T1 in the foil is equivalent to 26 /iBq/m 3 of 208 T1 in 
the gas volume. Figure [5] shows that this level of external 
background would lead toa~ 15% reduction in the half- 
life limit. The dominant 2z//3/3 background is measured by 
SuperNEMO and statistical uncertainties on its half-life 
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Fig. 4: Theoretical and experimental electron angular distributions for (a) MM and (b) RHCa- Theoretical and 
experimental electron energy difference distributions for (c) MM and (d) RHCa- All distributions are shown for the 
isotope 82 Se and the reconstructed distributions are normalised to the theoretical distribution to show signal efficiency. 



are expected to be negligible. Inclusion of an estimated 
7% correlated systematic uncertainty on the signal and 
background distributions [35] leads toa~ 5% reduction in 
the MM half-life limit. The effects of external background 
and of systematic uncertainties on the 2vf3(3 background 
are not included in the results of this paper. 

Expected exclusion limits at 90% confidence level on 
the half-life are shown in Fig.[7J Results are displayed as a 
function of RHCa admixture, where the signal distribution 
is produced by combining weighted combinations per bin 
of the MM and RHCa contributions at the event level. An 
admixture of 0% corresponds to a pure MM contribution, 
and an admixture of 100% to pure RHCa- Interference 
terms are assumed to be small and are neglected in the 
experimental simulation. The lower efficiency in the case 
of RHCa results in a lower limit for larger admixtures. 
The half-life limit is approximately twice as sensitive in 



measurements of 82 Se due to the lower mass number and 
higher 2v(3(3 decay half-life, though this is compensated 
in 150 Nd by more favourable phase space when convert- 
ing into physics parameter space. In the case where one 
mechanism dominates SuperNEMO is expected to be able 
to exclude 0v(5j3 half-lives up to 1.2 • 10 26 y (MM) and 
6.1 • 10 25 y (RHCa) for 82 Se, and 5.1 • 10 25 y (MM) and 
2.6 • 10 25 y (RHCa) for 150 Nd. 



4.3 Observation 

A Qv/3/3 signal rate with significant excess over the back- 
ground expectation, as for example shown in Fig. [5J would 
lead to an observation. The expected experimental sta- 
tistical uncertainties on the decay half-life are calculated 
from the Gaussian uncertainties on the observed number 
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Fig. 5: Expected number of MM signal (half-life of 10 25 y) and background events in 82 Se after 500 kg y exposure 
shown for (a) electron energy sum, (b) electron energy difference and (c) cosine of angle between electrons. 
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of signal and background events in the simulation. Fig- tistical uncertainty increases significantly for large admix- 

ure[5]shows the results for 82 Se and 150 Nd as a function turcs of RHCa at T x / 2 — 10 26 y which go beyond the 

of the admixture of RHCa- Acceptance effects cause the exclusion limit of SuperNEMO. 
uncertainty to increase with admixture of RHCa- The sta- 
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Fig. 8: One standard deviation statistical uncertainties in the measurement of double j3 decay half-lives at SuperNEMO 
as a function of admixture of the RHCa mechanism represented as band thickness for (a) 82 Se and (b) 150 Nd. 



The angular asymmetry parameter kg in Equation (|23p 
is experimentally accessible by defining N + as the number 
of events with measured angle cos 6 < and N_ as the 
number of events with cos 9 > 0. Similarly, an energy dif- 
ference asymmetry He can be obtained where N + is the 
number of events with energy difference < Q/2 (half the 
energy of the Oz//3/3 decay) and N_ is the number of events 
with energy difference > Q/2. The electron energy sum is 
required to be greater than 2.7 MeV for 82 Se and 3.1 MeV 
for 150 Nd to maximise signal to background ratio. This re- 
sults in signal efficiencies of 23.2% for the MM and 13.2% 
for the RHC A in 82 Se and 19.1% for the MM and 10.4% 
for the RHC A in 150 Nd. 

Experimentally, the distributions are only available as 
a sum of signal plus background so the measured values 
differ from the theoretically expected values due to the 
background distributions. This generally results in recon- 
structed correlation factors that are biased towards pos- 
itive values. The measured values of fcg,£ are shown in 
Fig. for a number of half-lives in the two isotopes. Sta- 
tistical uncertainties are shown as the width of the bands. 



All reconstructed kg y E values are displayed as a function of 
the corresponding theoretical kj E parameter, to allow for 
a model independent generalisation. It can be seen that 
the energy difference distribution allows stronger model 
discrimination than the angular distribution. 



5 Probing New Physics 

5.1 Model Parameter Constraints 

Having performed a detailed experimental analysis includ- 
ing a realistic simulation of the detector setup, the re- 
sults are interpreted in terms of the expected reach of the 
SuperNEMO experiment to new physics parameters of the 
combined MM and RHCa model of 0vf3(3 decay. 

Using Equation (fT3"]) for the Oz^/3/3 decay half-life to- 
gether with the coefficients listed in Table [TJ the expected 
90% CL limit on Tx/ 2 shown in Fig. [7] can be translated 
into a constraint on the model parameters m v and A. As- 
suming all other contributions are negligible this is shown 
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Fig. 9: Simulation of the correlation coefficients kg and kE as a function of theoretical kj E . The bands represent the 
one standard deviation statistical uncertainties. Shown are the angular correlation factor kg for 82 Se (a) and 150 Nd 
(b) and the energy difference correlation factor kE for 82 Se (c) and 150 Nd (d). 



in Fig. [TUJ (a) , as a contour in the m„ — A parameter plane. 
In the case SuperNEMO does not see a signal these pa- 
rameters would be constrained to be located within the 
coloured contour. The odd shape of the coloured contour 
is a direct consequence of the SuperNEMO 90% CL exclu- 
sion limit as a function of the specific admixture between 
the MM and the RHCa shown in Fig. [7] The small in- 
terference term, though not included in the experimental 
simulation, is taken into account through Equation (|13[) in 
this figure and results in the asymmetry of the distribution 
with respect to the sign of the parameter A. 

As shown in Section @J SuperNEMO is expected to 
be more sensitive to the Qvftfl half-life when using the 
isotope 82 Se, but this is compensated by the larger phase 
space of 150 Nd. As a result, the constraint on the model 
parameters is slightly stronger for 150 Nd. Due to the large 
uncertainty in the NMEs, this might be different for other 
NME calculations. To demonstrate the improvement over 
existing experimental bounds, the parameter constraints 
are shown in Fig.[TU](b) on a logarithmic scale (for positive 
values of A), comparing the expected SuperNEMO reach 



with the current constraints from the 0^/3/3 experiments 
NEMO-III [HI38] and Heidelberg Moscow [25]. 

Figure [TU] shows that SuperNEMO is expected to con- 
strain model parameters at 90% CL down to (m„) =70-73 
meV and A=(1-1.3)T0 . This would be an improvement 
by a factor 5-6 over the current best limit from the Heidel- 
berg Moscow experiment and more than an order of mag- 
nitude compared to the NEMO-III results. 



5.2 Angular and Energy Correlations 

As a more intriguing scenario it is now assumed that 
SuperNEMO actually observes a Oz//3/3 decay signal in 82 Se 
or 150 Nd. Because of the tracking abilities described in 
Section 2] this opens up the additional possibility of mea- 
suring the angular and energy distribution of the decays. 
Depending on the number of signal events detected, this 
can be crucial in distinguishing between different 0z//?/3 
decay mechanisms. In the analysis a reconstruction of the 
angular and energy correlation coefficients kg and kE is 
used to determine the theoretical coefficients, and thereby 
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the admixture between the left- and right-handed currents 
in the combined MM and RHCa model. 

For the isotope 82 Se, this is shown in Fig. [IT] for dif- 
ferent RHCa admixtures. The two blue elliptical contours 
correspond to the allowed one standard deviation (m v — A) 
parameter space at SuperNEMO when observing a signal 
at T 1/2 = 10 25 y and T 1/2 = 10 26 y, respectively. This 
takes into account a nominal theoretical uncertainty on 
the NME of 30% and a one standard deviation statistical 
uncertainty on the measurement determined from the sim- 
ulation (Fig. [5J. The blue elliptical error bands therefore 
give the allowed parameter region when only considering 
the total 0v/3(3 rate, which does not allow to distinguish 
between the MM and RHCa contributions. 

When taking into account the information provided 
by the angular and energy difference distribution shape, 
the parameter space can be constrained significantly. This 
is shown using the green contours in Fig. [TTJ for (a) a 
pure MM model, (b) a RHCa admixture of 30%, corre- 
sponding to an angular correlation of kg sa 0.4 and (c) 
a pure RHCa model. This fixes two specific directions 
in the m v — A plane (one for positive and one for neg- 
ative A). The widths of the contours are determined by 
the uncertainty in determining the theoretical correlation 
and admixture from the apparent distribution shape, see 
Fig. O The outer (light green) contours in Fig. [TTJ give 
the one standard deviation uncertainty on the parame- 
ters from reconstructing the angular distribution, while 
the inner (darker green) contour gives the one standard 



deviation uncertainty when using the distributions of the 
electron energy difference. As was outlined in Section [4j 
the energy difference distribution is expected to be easier 
to reconstruct and therefore gives a better determination 
of the RHCa admixture and a better constraint. While in- 
terference between MM and RHCa is neglected in the sim- 
ulation, it is taken into account in Equation f|13[) through 
the term C m \\i\ resulting in the slightly tilted elliptical 
contours and the asymmetry for A -H- — A. Finally, the red 
contours in Fig. [TTJ show the constraints on the model 
parameters when combining both the determination of 
the Ov/3/3 decay rate and the decay energy distribution. 
This demonstrates that such a successful combination can 
make it possible to determine the mechanism (i.e. the de- 
gree of MM and RHCa admixture in this case), and pro- 
vide a better constraint on the model parameters. From 
Fig. [TTJ (a), the Majorana mass term can be determined at 
(m v ) = 245+^ meV while the A parameter is constrained 
to be -0.87 • 10~ 7 < A < 0.92 • 10~ 7 in the case of a mea- 
sured 0vf3/3 decay half-life of 82 Se of T 1/2 = 10 25 y. For a 
82 Se half-life of T\i% = 10 26 y, the uncertainty on the decay 
rate increases as SuperNEMO reaches its exclusion limit 
for RHCa admixtures. It is therefore only possible to ex- 
tract upper limits on the model parameters from Fig. [TTJ 
for X1/2 = 10 26 y. However, the shape information pro- 
vides additional constraints on the parameter space. In 
Fig. [T^] we show the analogous plots for the isotope 150 Nd 
assuming a decay half-life of Ti/ 2 = 10 25 y. 
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Fig. 11: Constraints at one standard deviation on the model parameters m v and A for Se from: (1) an observation 
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300 
250 

>200 

S 
X150 

v ioo 

50 





300 
250 

>200 

E 

X150 
S 
^100 

50 






300 




250 




>200 

E 

"150 

> 




100 




50 


~~~~~~~i 





1 




Fig. 12: As Fig. [TT] but for the isotope i5U Nd with a decay half-life of T 1/2 = 10 25 y 



5.3 Rate Comparison of 150 Nd and 82 Se 

While reconstruction of the decay distribution can be an 
ideal way to distinguish between different mechanisms, it 
might be of little help if Oi//3/3 decay is observed close to 
the exclusion limit of SuperNEMO, or not at all. This is 
demonstrated in Fig. [TT] where, for a half-life of T x / 2 — 
10 26 y, the reconstruction of the energy difference distri- 
bution will be problematic due to the low number of events 
(compare Fig.[3]). As an alternative, it is possible to com- 
pare the 0vj3f3 rate in different isotopes. This method, 
which could provide crucial information close to the ex- 
clusion limit, is especially relevant for SuperNEMO which 
could potentially measure 0v/3(3 decay in two (or more) 



isotopes. Such a comparative analysis was used in |21j 
to distinguish between several new physics mechanisms. 
A combined analysis of several isotopes, potentially mea- 
sured in other experiments, will improve the statistical 
significance |22) . 

The possibility of sharing the two isotopes equally in 
SuperNEMO, each with a total exposure of 250 kg y, is 
now considered. In the cases where the MM or the RHCa 
contributions dominate, the following half-life ratios can 
be found: 



MM 



T s2 Sc 
J_L/2_ 

T 150 Nd 

■n/2 



a 



Nd 
mm 



(2-7) 2 'C 



2.45, 



(27) 
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T 82 Sc 

J l/2 
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£< 15(1 Nd 



(2.7) 



2 r> 82 Se 



3.64. 



(28) 



These ratios and their uncertainties are determined by the 
0vj3f3 decay NMEs and phase spaces. The factor 2.7 is the 
correction added to the 150 Nd NMEs as described in Sec- 
tion 12.31 It has recently been suggested that uncertainties 
in NME calculations are highly correlated [43 so mea- 
surements made with two or more isotopes could reduce 
the uncertainty on the physics parameters significantly. 
Additionally, most experimental systematic uncertainties 
would cancel if different isotopes are studied in a single 
experiment such as SuperNEMO. This would not be pos- 
sible when comparing results with other experiments. The 
statistical uncertainties are naturally greater than in the 
single-isotope case, due to the exposure being halved for 
each isotope. 

The results of the combined NME and statistical un- 
certainties analysis, including a possible correlation of the 
NMEs, are illustrated in Fig. [£31 It shows the Oiy/3/3 half- 
life of 150 Nd as a function of the half- life in 82 Se assuming 
a pure MM model, with the coloured contours giving the 
deviation from the hypothesis that the mass mechanism 
is the single source of 0vj3/3 decay in both isotopes at the 
1, 2 and 5 standard deviation level. The statistical uncer- 
tainties used in Fig. ll3l are derived from our experimental 
simulation and the standard 30% NME uncertainties are 
applied. The effect of a possible correlation of the NMEs 
is shown by assuming the covariance coefficient between 
the NME uncertainties of 82 Se and 150 Nd to be (a) zero 
(no correlation), (b) 0.7 and (c) 1.0 (full correlation). The 
experimental uncertainties and expected sensitivity (90% 
CL exclusion) limits are calculated for 250 kg y of ex- 
posure of each isotope and assume a 50% 82 Se and 50% 
150 Nd option for SuperNEMO. The red line shows the re- 
lationship for the half-life ratio in the pure RHCa model 



(Equation (|2"8")) ). It can be seen that an exclusion at two 
standard deviations is possible if the NME errors are per- 
fectly correlated and at the one standard deviation level if 
the correlation is 70%, which is a more realistic assump- 
tion. 

Other mechanisms have different half-life ratios [2T] 
so they could be excluded with different CLs at Super- 
NEMO. One important advantage of this method is that 
it provides a possibility to falsify the mass mechanism as 
the sole source for 0vfif3. A measurement within the blue 
contour would indicate that the pure MM model can be 
excluded at the 5 standard deviation level and new physics 
is required to explain the measured half-lives. 



5.4 Combined Energy and Rate Comparison of 150 Nd 
and 82 Se 

In the most favourable case, signal event rates in two iso- 
topes could be high enough (0^/3/3 decay half-lives small 
enough) that the distribution method and the two iso- 
tope rate analysis can be combined to put further con- 
straints on the parameter space. The effect of such a com- 
bined analysis on the allowed parameter space is shown 



in Fig. [H where the 50% 150 Nd - 50% 82 Se two-isotope 
option (red contours) is compared to the single-isotope op- 
tions 100% 82 Se (green contours) and 100% 150 Nd (blue 
contours). The 0^/3/3 decay half- life of 82 Se is assumed 
to be 10 25 y and the half-life of 150 Nd is determined 
by the respective MM-RHC A admixture, i.e. (a) T^ = 

10 25 /2.45 y, (b) 10 25 /2.73 y and (c) 10 25 /3.64 y. The NME 
uncertainties are assumed to be 30% with a 0.7 covariance 
between the uncertainties of the NMEs of 82 Se and 150 Nd. 
As can be seen in Fig. 1141 the two-isotope option can im- 
prove the constraints on the parameter space along the 
radial direction, e.g. it allows a more accurate determina- 
tion of the MM neutrino mass m v in Fig. [14] (a). On the 
other hand, the accuracy in the lateral direction (the pa- 
rameter A in Fig. [2](a)) becomes worse compared to the 
best single-isotope option due to the reduced statistics for 
a given isotope. 



6 Conclusion 

The 0^/3/3 decay is a crucial process for physics beyond the 
Standard Model, and the next generation SuperNEMO ex- 
periment is designed to be a sensitive probe of this lepton 
number violating observable. In addition to being able to 
measure the 0i//3j3 half-life of one or more isotopes, it also 
allows the determination of the angular and energy differ- 
ence distributions of the outgoing electrons. 

In this paper we have focussed on the sensitivity of 
SuperNEMO to new physics and its ability to discriminate 
between different §v(3(3 mechanisms. This was achieved by 
a detailed analysis of two important models, namely the 
standard mass mechanism via light left-handed Majorana 
neutrino exchange and a contribution from right-handed 
current via the effective A parameter stemming from Left- 
Right symmetry. The study included a full simulation of 
the process and the SuperNEMO detector at the event 
level, allowing a realistic estimation of the experimental 
90% CL exclusion limit and statistical uncertainties. 

SuperNEMO is expected to exclude 0i//3/3 half-lives up 
to 1.2 • 10 26 y (MM) and 6.1 • 10 25 y (RHC A ) for 82 Se and 
5.1 ■ 10 25 y (MM) and 2.6 • 10 25 y (RHC A ) for 150 Nd at 90% 
CL with a detector exposure of 500 kg y. This corresponds 
to a Majorana neutrino mass of m„ s=s 70 meV and a A 
parameter of A w 10~ 7 , giving an improvement of more 
than one order of magnitude compared to the NEMO-III 
experiment. 

It has been shown that the angular and electron en- 
ergy difference distributions can be used to discriminate 
new physics scenarios. In the framework of the two mecha- 
nisms analysed, it was demonstrated that using this tech- 
nique the individual new physics model parameters can be 
determined. For a half-life of T x / 2 = 10 25 y with an expo- 
sure of 500 kg y, the Majorana neutrino mass can be de- 
termined to be 245 meV with an uncertainty of 30% while 
the A parameter can be constrained at the same time to 
be smaller than |A| < 0.9 • 10~ 7 . Such a decay distribution 
analysis could be easily extended further to include other 
new physics scenarios with distinct distributions and the 
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Fig. 13: The 0vf3(3 half-life of 150 Nd as a function of measured half-life in 82 Se for the hypothesis that the MM 
is the single Ov/3/3 decay source. The contours show the 1, 2 and 5 standard deviation levels assuming statistical 
uncertainties derived from the experimental simulation and 30% NME errors assumed to have (a) no, (b) 0.7 and (c) 
perfect correlation. The experimental uncertainties and expected sensitivity (90% CL exclusion) limit are calculated 
for 250 kg y of exposure (assuming a 50% 82 Se and 50% 150 Nd option). The red line shows the relationship for the 
RHCa- The blue contour shows the 5<r exclusion of the MM. 
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Fig. 14: Constraints at one standard deviation on the model parameters m„ and A from: (1) an observation of Of/3/J 
decay half-life of 82 Se at l\/ 2 = 10 25 y with 500 kg y exposure and reconstruction of the energy difference distribution 
(outer green contour); (2) an observation of 0v/3j3 decay half- life of 150 Nd at a half- life corresponding to 2\ « = 10 25 y 
in 82 Se with an exposure of 500 kg y and reconstruction of the energy difference distribution (inner blue contour); 
(3) combined analysis of (1) and (2) with an exposure of 250 kg y in 82 Se and 150 Nd (red contour). The admixture 
of the MM and RHCa contributions is assumed to be: (a) pure MM contribution; (b) 30% RHCa admixture; and (c) 
pure RHCa contribution. NME uncertainties are assumed to be 30% with a correlation of the uncertainties of 0.7, and 
experimental statistical uncertainties are determined from the simulation. 



results are quoted in terms of a generalised distribution 
asymmetry parameter to allow new physics scenarios to 
be compared. As the two example mechanisms considered 
exhibit maximally different angular and energy distribu- 
tion shapes, they serve as representative scenarios cover- 
ing a large spectrum of the model space. For example, 
the right-handed current contribution due the effective rj 
parameter, also arising in Left-Right symmetrical models, 



can be distinguished from the mass mechanism and the 
right-handed current A contribution by looking at both 
the angular and energy difference decay distribution. This 
would allow a determination of all three model parameters 
m v , A and rj by looking at the total rate and the angular 
and energy difference distribution shapes. 

Further insight into the mechanism of 0^/3/3 can be 
gained by using multiple isotopes within the SuperNEMO 
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setup. This possibility was explored by studying the op- 
tion of having 50% 150 Nd and 50% 82 Se, each with an 
exposure of 250 kg y While the statistics per isotope is re- 
duced compared to the individual 100% options, the abil- 
ity to measure the ratio between the half-lives of the two 
isotopes can be used as additional information on the un- 
derlying physics mechanism responsible for Oz^/3/3 decay. 
As was shown for the isotopes 82 Se and 150 Nd at Super- 
NEMO, this could be a powerful method to falsify the 
mass mechanism as the dominant 0^/3/3 mechanism. A cor- 
relation between the uncertainties of nuclear matrix ele- 
ments, which is generally expected on theoretical grounds, 
has proven to be of importance and its impact on the falsi- 
fication potential was analysed. Within SuperNEMO such 
a correlation could also be found between the systematic 
uncertainties in the measurements of different isotopes. 

SuperNEMO also has a number of other possibilities 
to disentangle the underlying physics. The detection tech- 
nology is not dependent on one particular isotope and any 
double j3 emitting source could be studied in the detector. 
In this paper 82 Se and 150 Nd have been considered but 
other isotopes such as 48 Ca or 100 Mo are feasible. The 
analysis can be extended to cover more than two isotopes 
thereby achieving a higher significance and a comparison 
with other experimental results will provide additional in- 
formation. SuperNEMO is also able to measure a Oz^/3/3 
decay to an excited state, by measuring two electrons and 
an accompanying photon. This again could be used to aid 
the analysis to discriminate between new physics mecha- 
nisms. 

A combination of the above methods makes Super- 
NEMO an exciting test of new physics. These methods 
would prove invaluable in excluding or confirming domi- 
nating mechanisms of lepton number violation in the reach 
of the next generation 0^/3/3 experiments. 
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